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Abstract 

An alternative representation of an exact statistical field theory for simple fluids, 
based on the method of collective variables, is presented. The results obtained are 
examined from the point of another version of theory that was developed recently 
by performing a Hubbard-Stratonovich transformation of the configurational Boltz- 
mann factor [J.-M. Caillol, Mol. Phys. 101 (2003) 1617]. The analytical expressions 
for the pressure and the free energy are derived in two-loop approximation for both 
versions of theory and it is shown that they are indeed equivalent. The results yield 
a new type approximation within an untested approximation scheme. 
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1 Introduction 



Functional methods in modern statistical physics represent one of the most 
powerful tools for the study both of equilibrium and dynamical properties 
(see, e.g. [1,2]). A great amount of statistical field theories known in the liter- 
ature are based of the Hubbard-Stratonovich transformation [3,4], proposed in 
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the 50ies. Nearly at the same time another method - so-called collective vari- 
ables (CV) method - that allows in a explicit way to construct a functional 
representation for many-particle interacting systems was developed [5,6] and 
applied for the description of charged particle systems, in particular, to the 
calculation of the configurational integral of the Coulomb systems. The idea 
of this method was based on: (i) the concept of collective coordinates being 
appropriate for the physics of system considered (see, for instance, [7]), and 
(ii) the integral equality allowing to derive an exact functional representation 
for the configurational Boltzmann factor. Later the CV methods was success- 
fully developed for the description of classical many-particle systems [8] and 
for the phase transition theory [9]. One more functional approach, the meso- 
scopic field theory, was recently developed for the study of phase transitions 
in ionic fluids [10]. 

One of the goals of this paper is to reconsider the CV method from the point 
of view of the statistical field theory and to compare the results obtained 
with those found recently by one of us by means of the KSSHE (Kac-Siegert- 
Stratonovich-Hubbard-Edwards) theory [11]. 

We formulate the method of CV in real space and consider a one-component 
continuous model consisting of hard spheres interacting through additive pair 
potentials. The expression for the functional of the grand partition function 
is derived and the CV action that depends upon two scalar fields - field p 
connected to the number density of particles and field u conjugate to p - is 
calculated. We study correlations between these fields as well as their relations 
to the density and energy correlations of the fluid. The grand partition func- 
tion of the model is evaluated in a systematic way using a well-known method 
of statistical field theory, namely the so-called loop expansion. It consists in 
expanding functionally the action H around a saddle point, so that the low- 
est order (zero loop) approximation defines the mean-field (MF) level of the 
theory and the first order loop expressions correspond to the random phase 
approximation (RPA). Recently [11] this technique was applied to the action 
obtained within the framework of the KSSHE theory. In this paper we perform 
a two-loop expansion of the pressure and the free energy of the homogeneous 
fluid which yields a new type of approximation which we plan to test in our 
future work. 

The paper is organized as follows. In Section 2, starting from the Hamiltonian, 
we introduce the two different functional representations of the grand partition 
function based on the KSSHE and CV methods. Here we also enter several 
types of statistical field averages that are important in the further part of 
the paper. In Section 3 we introduce the CV and KSSHE field correlation 
functions, establish links between them as well as their relation to the density 
correlation functions of the fluid. The MF level of the KSSHE and CV field 
theories is formulated in Section 4. Section 5 is devoted to the loop expansion 
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of the grand potential. The pressure and the free energy of the homogeneous 
fluid are obtained in the two-loop approximation in Section 6. We conclude 
with some final remarks in Section 7. 



2 Summary of previous works 

2.1 The model 

We consider the case of a simple three dimensional fluid that consists of iden- 
tical hard spheres of diameter a with additional isotropic pair interactions 
v(rij) [r\j = |xj — Xj|, Xj is the position of particle "i"). Since v(r) is arbi- 
trary in the core, i.e. for r < a, we assume that v(r) has been regularized 
in such a way that its Fourier transform v q is a well behaved function of q 
and that v(0) is a finite quantity. We denote by Q the domain of volume V 
occupied by the molecules of the fluid. The fluid is at equilibrium in the grand 
canonical (GC) ensemble and we denote by (3 = 1/kT the inverse tempera- 
ture (k is the Boltzmann constant) and /x is the chemical potential. In addition 
the particles are subject to an external potential ip(x) and we will denote by 
f(x) = /3(/a — ipi'x)) the dimensionless local chemical potential. We will stick 
to notations usually adopted in standard textbooks on the theory of liquids 
(see e.g. [12]) and denote by w(r) = —f3v{r) the negative of the dimensionless 
pair interaction. Quite arbitrarily we will say that the interaction is attractive 
if the Fourier transform w(q) is positive for all q; in the converse case it will 
be said repulsive. 

In a given configuration C = (N; Xi . . . xjv) the microscopic density of particles 
reads 

N 

p(x|C) = E^ (3) (x-x l ) ) (1) 
1=1 

and the GC partition function E [u] can thus be written as 

00 1 r f 1 \ 

s M = E L dl • • ■ dn 6XP {-^^ + 2 {Plwm + W) ) ' (2) 

where % = Xj and di = d 3 Xi. For a given volume V, S [i/] is a function of /3 
and a convex functional of the local chemical potential v(x) which we have 
strengthened by using a bracket. In eq. (2) exp(— PVus(C)) denotes the hard 
sphere contribution to the Boltzmann factor in a configuration C and V = u+u$ 
where z/g = — w(0)/2 is proportional to the self-energy of the particles. From 
our hypothesis on w(r), v$ is a finite quantity which depends however on the 
regularization of the potential in the core. In the r.h.s of eq. (2) we have also 
introduced Dirac's brac-kets notations 
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(p\ 



(v\p) = fdi nmi) 

Jn 

\w\p) = [ dld2 p(l\C)w(l,2)p(2\C) . 



Jn 



(3b) 



(3a) 



Previous work have shown that E [is] can be rewritten as a functional integral. 
Two such equivalent representations are reviewed below. 

2.2 The KSSHE representation 

As it is well known the GC partition function E [u] can be re-expressed as a 
functional integral by performing the KSSHE transformation [13,14,4,3,15,16] 
of the Boltzmann's factor. Under this transformation E [u] can be rewritten 
as the GC partition function of a fluid of bare hard spheres in the presence 
of a random Gaussian field tp with a covariance given by the pair potential 
[17,19,11,18]. More precisely one has 

• i) in the attractive case (w(q) > 0) 



where, in both cases, <p is a real random field and Shs denotes the GC partition 
function of bare hard spheres. E can thus be written as a Gaussian average 
(. . .) w of covariance w and we have noted by J\f w the normalization constant 



The functional integrals which enter eqs. (4), (5) and (6) can be given a precise 
meaning in the case where the domain Q is a cube of side L with periodic 
boundary conditions (PBC) which will be implicitly assumed henceforth. More 
details are given in Appendix A. 

In the repulsive case, the hard core part of the interaction is not compulsory 
for the existence of a thermodynamic limit [20] and the reference system can 
be chosen as the ideal gas [21,22]. 

Eqs (4) and (5) are easily generalized to the case of molecular fluids or mix- 
tures, for instance a charged hard spheres mixture (primitive model) [18,23]. 




(4) 



E [v] =JV { _! W) j Vip exp ^- (tp\w V)j s hs [v + M 



(5) 
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When the pair interaction w is neither attractive nor repulsive it is necessary 
to introduce two real scalar fields if some rigor is aimed at [11]. Alternatively, 
eq. (4) can be considered to hold in any case having in mind that <p can be a 
complex scalar field. Therefore we shall write formally in all cases 

2 M = K 1 J Vip exp (-W K [v, <p}) , (7) 
where the action of the KSSHE field theory reads as 

H K W,<p] = ^ (</?|w~V) - lnS HS [v+ <p] . (8) 
2.3 The CV representation 

We introduce now briefly the CV representation of S [v] and refer the reader 
to the literature for a more detailed presentation [9,6,8,24,25]. The starting 
point is the identity 

exp Q (pMp)) = J Vp5 r [p- p] exp Q (p|Hp)) 

= | VpVu exp Q (pMp> + i H {p - p})) , (9) 

where we have made use of the functional "delta function" [1] 5? [p] defined 
in eq. (A. 6) in appendix A. Inserting eq. (9) in the expression (2) of the GC 
partition function one finds 

E [z/] = J VpVu exp (-Wcv [u, p, u}) , (10) 

where the action of the CV field theory reads as 

Hcv[v,P,u] = ~{p\w\p) -i(u\p) -lnS HS [^-H . (11) 

We stress that u and p are two real scalar fields and that eqs. (10) and (11) are 
valid for repulsive, attractive as well as arbitrary pair interactions. Moreover, 
with the clever normalization of Wegner [19] for the functional measures there 
are no unspecified multiplicative constant involved in eq. (10) (see Appendix 
A for more details). 

The CV transformation is clearly more general than the KSSHE transforma- 
tion since it can be used for a pair interaction w(l, 2) which does not possess 
an inverse and is easily generalized for n-body interactions with n > 2. The 
equivalence of the CV and KSSHE representations (7) and (10) of S [v] is 
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readily established in the repulsive case ( w < 0) by making use of the prop- 
erties of Gaussian integrals (cf. eq. (A. 5) of Appendix A). In the attractive or 
in the general case we did not find a convincing way to obtain one formula 
from the other. 



2.4 Statistical average 

In the sequel it will be important to distinguish carefully, besides the usual 
GC average < A(C) >gc of a dynamic variable A(C), between several types 
of statistical field averages. Firstly Gaussian averages of covariance w 

{A[<P]) v ,=K 1 fv<P exp(-i(v?KV))-4M, (12) 

where A [ip] is some functional of the KSSHE field <p and M w has been defined 
in eq. (6); secondly the KSSHE averages defined as 

(A M) K = 2 W 1 Jv^A [if] exp (-Hk [u, <p\) , (13) 

and thirdly the CV averages defined in a similar way as 

(A [p, lu]) cv = E [vY 1 j VpVLU A [p, u] exp (-Hcv [v, P, v]) , (14) 

where A [p, uj] is a functional of the two CV fields p an uj. 



3 Correlation functions 



Since all thermodynamic quantities of the fluid can be expressed in terms of 
the GC density correlation functions [12] it is important to relate them 
to the field correlation functions in the CV and the KSSHE representations; 
this is subject of section 3. 



3.1 Density correlations 

The ordinary and truncated (or connected) density correlation functions of 
the fluid will be defined in this paper as [12,26,27] 
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\i=i / GC S M 6v(l)...6v(n) ' 

Our notation emphasizes the fact that the (connected and not connected) 
are functional of the local chemical potential v(x) and functions of the coor- 
dinates (1, . . . , n) = (xi, . . . , x n ). 

We know from the theory of liquids that [26,27] 

G lB) ' T M(l. ■■■,") = G (B) H(1, ..,n)-EII G (m) ' T M(^i, • • • ,<m) , (16) 

where the sum of products is carried out over all possible partitions of the set 
(1, . . . ,n) into subsets of cardinal m < n. Of course p[z/](a;) = G^ n=1 ^[u](x) = 
G( n=1 ^ T [v\(x) is the local density of the fluid. 

It follows from the definition (15) of the G^[z/](1, . . . ,n) that they can be 
reexpressed as KSSHE or CV statistical averages, i.e. 



G^[u}(l,...,n) = (G^[v+ V ](l,...,n)) K , (17a) 
G^[u](l,...,n) = ((^[v-iu](l,...,n)) cv . (17b) 

Although enlightening these relations are not very useful except for the special 
case n = 1 which reads explicitly as 



p[v](x) = {pHs[v + <p](x)) K , (18a) 
PW) ( x ) = (Phs[^-M( x ))cv > ( i8b ) 

where pns[£]( x ) is the local density of the hard sphere fluid with the local 
chemical potential £(x). 



3.2 KSSHE field correlations 



Let us introduce the modified partition function 

- 1 [u, J] = AT 1 J TV exp (-H K [v, ¥>] + (A<P)) , (19) 

where J is a real scalar field. Clearly S 1 [u, J] is the generator of field correla- 
tion functions and a standard result of statistical field theory yields [2] 
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Moreover one has [2] 



j=o 

(20) 



j=o 



(?W' T (1, . . . , n) = G*">(1, . . . , n) - £ II G^' T (ii, . . . , v») . (21) 



m<n 



The relations between the G^' T and the truncated density correlation func- 
tions G^' T have been established elsewhere by one of us [11,18,23]. We sum- 
marize them below for future reference. 



pW](l)=w-\l,l)(v(l)) K , (22a) 
G™' T [v] (1, 2) = -w-\l, 2) + w-\l, 1>~ 1 (2, 2) x 

xG? r H(l',2), (22b) 
G in) ' T [u] (1, . . . , n) = w-\l, I)... w-\n, n) x 

x G^ ] ' T [u] (1', . . . , n') forn > 3 , (22c) 

where we have adopted Einstein's convention, i.e. a space integration of posi- 
tion variables labled by the same dummy indices over the domain Q is meant. 

3.3 CV field correlations 



In this section we study the correlation functions of the fields p and uj in the 
CV representation. We thus define 



Gf)[i/](l ) ..,n)= n^x*) , 
\i=i / cv 

G^M(l,...,n) = (n^(x i )\ , (23) 
\i=i / cv 

and their connected parts 



GW- r (l, . . . , n) = G?)(l, , n) - £ II >' r (ii, . . . , i m ) , 

Gi n) ' T (l, . . • , n) = Gf)(l, ..,n)-EIl G^ft, . . . , U . (24) 



m<n 
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3. 3. 1 Correlation functions 

Let us define the modified GC partition function 



E 2 [u, J] = J VpVu exp (-Wcv [v, P, + W) 



(25) 



where J is a real scalar field. H 2 [z/, J] is clearly the generator of the and 
we thus have 



^ W'-"'^ S [i/] 5J(1) . . . 5 J(n) 

c (n) ' T M(i n) = i^!!^L 

^ ^'•••'^ 5J(l)...5J(n) 



j=o 



(26) 



j=o 



The simplest way to obtain the relations between the G^ and the density 
correlation functions is to start from the definition (15) of G^ n \ One has 



G (n)(l n) = J S M 



^jVp-Du exp G<HHp)+^W) ^y^ 

^ | DpDo; exp Q + ^ (o;|p)) x 

n <5 n Shs[^-H 



X W (Jw(l)...M«) ' 
Performing now n integral by parts yields 



G (n) (l,...,n) = m^ J VpVu exp Q (p\w\p) + In S HS [^ - ^]) x 
5™ exp (i (o;|p)) 



U=l / cv 



5u;(l) . . . 5a; (n) 
We have thus proved the expected result 

G (n) H (l,...,n) = (?W [i/](l,...,n), 
which implies in turn that 

G^> T [v] (l,...,») = Gf M(l,...,n). 



(27) 



(28) 
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3.3.2 Correlation functions G^> 

Let us define the modified GC partition function 

~ 3 [u, J] — J VpVu exp (-Hcv W, P, to] + (J\u)) , 



(29) 



where J is a real scalar field. H 3 [v, J] is the generator of the and we thus 
have 



(#>[i/](l,...,n) = - 



5 n Z s [v, J] 



GW- r [i/](l,...,n) 



H[i/] <?J(1) . . . <SJ(n) 
<5 n log S 3 [i/, J] 



j=o 



<JJ(l)...5J(n) 



(30) 



j=o 



In order to relate Gffl to G^ we perform the change of variable p — > p + i J 
in eq. (29). The functional Jacobian of the transformation is of course equal 
to unity and one obtains the relation 



lnS 3 [u, J] = -- (J\w\J) + lnS 2 [v,iw*J] , 



(31) 



where the star * means a space convolution and H 2 was defined at eq. (25). 
The idea is to perform now n successive functional derivatives of both sides of 
eq. (31) with respect to J. Since it follows from the expression (26) of G^' T 
that 



S n log S 2 [is, iw * J] 



SJ(l)...SJ(n) 
one readily obtains 



j=o 



i n w(l, l') . . .w(n,n) x 
xG^ T [u](l',...,n), 



(32) 



Ml)) cv = ^(l,l')(p(l')) cv , (33a) 

Gg)' r [i/] (1,2) = -w(l, 2) 

- W (l,l>(2,2')Gf>](l',2'), (33b) 

G£ n) ' T [is] (1, . . . , n) = i n l') . . . iu(ra, n') x 

x G^' T [i/] (1', . . . , n') for n > 3 . (33c) 

A comparison between eqs (22), (28), and (33) will show that the CV and 
KSSHE correlation functions are simply related since we have, as expected 



G^ T H (1, . . . , n) = (-i) n G^> T [u] (1, . . . , n) for all n . (34) 
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4 Mean Field theory 



4-1 KSS HE representation 

We summarize previous work on the mean-field (MF) or saddle point approx- 
imation of the KSSHE theory [11,18]. At the MF level one has [2] 

2 M f H = exp (-H K [v, <A)]) , (35) 
where, for ip = <p , the action is stationary, i.e. 

5 H K [v, if] 



5(p 



= . (36) 



<P0 



Replacing the KSSHE action by its expression (8) in eq. (36) leads to an 
implicit equation for ip : 

V?o(l) = ^l,l')pHs[^ + y?o](l'), (37) 

which reduces to 

ifo = w(0) pns P + (Po] (38) 
for a homogeneous system. It follows from the stationary condition (36) that 
the MF density is given by 

Pmf M (1) = ^ n ^ M = Phs [V + <p ] (1) , (39) 
and that the MF grand potential reads 

lnS MF H = lnS HS p+ Vol ~ ^ (Pmf\w\pmf) ■ (40) 
Moreover, the MF Kohn-Scham free energy defined as the Legendre transform 
(3Amf [p] = sup {(p\u) - In S MF M} (41) 

is found to be 

PAmf [p] = PAhs [p] - \ (p\w\p) + l/ n dx w(0)p(x) . (42) 

It can be shown [11] that ^4mf [p] constitutes a rigorous upper bound for the 
exact free energy A [p] if the interaction is attractive (w(q) > 0) and a lower 
bound in the converse case (w(q) < 0). 

Finally, the pair correlation and vertex functions at the zero-loop order which 
are defined respectively as 
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G - H(i,2) - M1)M2) , 

n {2) r, 9 x 8 2 ft 'AMF [P] 

C MF[p}(^) = - 5pil)5pi2 y 



5 2 In S MF \v 



are given by 



(43a) 
(43b) 



G 



2 r 'l, 2) = (l-w G§ T [17 + ^ ■ G ( if [V + ipo] (1, 2) 



MF 



cS(l, 2) = -Ggf -'(1, 2) = C HS (1, 2) + w(l, 2) . 
It follows then from eq. (22b) that we have 

G^Sf(1, 2) = (l-w G ( il \V + <p ]) ^ ■ w(l, 2) 



(44a) 
(44b) 

(45) 



4-2 CV representation 



An analysis similar to that of Sec 4.1 can be made in the CV representation. 
The MF level of the CV field theory will be defined by 



2 M f H = exp (-Wcv [v, Po, wo]) , 
where, for p = p an d ^ = the CV action is stationary, i.e. 



5 ?icv [v, P, <A 



5p 



5 H C v [v, p, u] 



Sou 



. 



(46) 



(47) 



(POM)) 



Replacing the CV action by its expression (11) in eq. (47) yields a set of two 
coupled implicit equations for p and uq: 



= W (l,2)p (2)+^o(l) , 
= p (l) - Phs p-iuo] (1) • 



(48) 



If we define (p = —iujq, then the two previous equations can be rewritten 
ip = w(l,2)p (2) , 

Po(1)=Phs[^ + <A)](1) , (49) 

which shows that, as expected, ipo coincides with the saddle point of the 
KSSHE field theory (cf Sec 4.1). Moreover a direct calculation will show that 



In S M f M = -Hk [v, <Po] = ~^cv [v, Po, wo] 



(50) 
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Therefore the local density, the grand potential, the free energy, the correlation 
and vertex functions coincide at the MF level in the CV and KSSHE field 
theories. 



5 Loop expansion of the grand potential 

A one loop expression for the grand potential In S [u] and the free energy f3A [p] 
has been obtained in [11] in the framework of KSSHE field theory. Here we 
perform a two-loop expansion for InS. The derivation will be made both in 
the KSSHE and CV representations. 



5. 1 KSSHE representation 

The loop expansion of the logarithm of the partition function of a scalar field 
theory can be found in any standard text book, see e.g. that of Zinn- Justin 
[2]. One proceeds as follows. A small dimensionless parameter A is introduced 
and the loop expansion is obtained by the set of transformations 

y? = y?o + A 1/2 x , 

= In S<°> [u] + A In H« [u] + A 2 In [u] + C(A 3 ) , (51) 

where (f is the saddle point of the KSSHE action. At the end of the calculation 
one set A = 1. It follows from the stationary condition (36) that a functional 
Taylor expansion of Hk \y, V 9 ] about the saddle point has the form 

[v, f\ = H K [u, y? ] 1 / . A -i| \ 
A A 2 \ Xl volX / 

+ / dl...dnWW(l,...,n) X (l)...x(n), (52) 

n=3 U - JQ 

where we have adopted Zinn- Justin's notations [2]. In eq. (52) 

A-^l, 2) = w~\l, 2) - Gg- T [V + v^o] (1, 2) • (53) 

A vo (l, 2) is the free propagator of the theory, it coincides with C^ MF (1, 2) as 
can be seen by comparing eqs. (45) and (53). If v is uniform then the system 
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is homogeneous and A v , (l, 2) takes on a simple form in Fourier space, i.e. 



w{q) 

1 - ti(q)Ggl' J > + tpo] (q) 



K* (?) = ; ^wp- — • (54) 



Finally the kernels Tif$ in the r.h.s. of eq. (52) are given by 



H% (1, . . . , n) = -G<g« T [17 + v^o] (1, • • • , n) . (55) 



The expression of In [v] in terms of the propagator A 1 (1, 2) and the vertex 
interactions TC^ is obtained by means of a cumulant expansion of In S [v\ and 
by making a repeated use of Wick's theorem. 

Of course [i/] = Smf M- At the one-loop order one finds [11] 

=<■) H = ^ = J^^HHA>), m 

N w f Dtp exp (—1 (ip\w\(p)j 

For a homogeneous system the Gaussian integrals in the r.h.s. of eq. (56) can 
be performed explicitly (cf. eqs. (A. 5) of Appendix A) and one has 

In SW [u] = ~ jf In (l - w(q)G%l> T [V + <po] (?)) • (57) 

As is shown in detail in [11] the one-loop approximation coincides with the 
well known RPA approximation of the theory of liquids [12]. 

The second-loop order contribution lnS( 2 ) to the grand potential has a com- 
plicated expression involving the sum of three diagrams sketched in fig. 1 

\nZW[v] = D a + D b + D c . (58) 




Da D b D c 



Fig. 1. Diagrams which contribute to lnH^ 2 ) [v\. D a and D c are irreducible while 
Db is reducible. 



More explicitly one has [2] 
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D a = - [ rfl...rf4A w (l,2)A w (3,4)GS T [l7+^o](l,2,3,4) 

O JQ 

D h = \ [ dl...d3dl ...d3 A W) (l,2)A W) (l',2')A W) (3,3') 

o JQ 

x T [17 + ^o] (1, 2, 3)G§ T [17 + (1', 2', 3') 
D c = —f dl...d3dl' ...d3' A w (l,l')A w (2,2')A V0 (3,3') 

x T [V + (1, 2, 3)G@' T [17 + ^o] (1', 2', 3) . (59) 

As they stand, the above relations are not particularly useful for practical 
applications (even in the homogeneous case) since they involve the three and 
four body correlation functions of the reference HS fluid. We will introduce 
some reasonable approximation in Sec. 6 to tackle with this horible expression. 
Quite remarkably it has been shown recently that for a symmetric mixture of 
charged hard spheres In has a much more simple expression which involves 
only the pair correlation functions of the HS fluid as a consequence of local 
charge neutrality [18,23]. 



5.2 CV representation 



In order to obtain a loop expansion of In S in the CV representation we con- 
sider the following set of transformations 



p = p + X 1/2 Sp , 



In S [v] — A In i^J V5p V5uo exp 



Hcv [v, p, to] 
X 



(60a) 
(60b) 

(60c) 



where (po,^o) is the saddle point of the CV action. The form retained in 
eqs (60a) and (60b) is imposed by the exact relation (uj) cv = iw * (p) cv 
derived in Sec 3.3.2. It follows from the stationary condition (47) that the 
functional Taylor expansion of TLc\ [v, p, u] about the saddle point reads as 



Hcv[v,P,w] ^cv[^,Po,^o] 1 /x i ,<- v -,x \x\ 
^ = ^ 2 ( S P\ W \ 6 P) ~ 1 \ 5uj \ 5 P) 



1 



5uj 



G (2), T Ft7 



6u) + 6H[6u] , 



(61) 



r H g [1/-IU \ 

where 

5H[5lu] = -J2 ^J-X^" 1 J n dl...dnH^(l,...,n)5uj(l)...5uj(n) , (62) 



n=3 
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since, as </?o = — i^o we have indeed Cr^g T \ v ~ ^o] = where the kernel 

is precisely that defined in eq. (55). 

We are thus led to define a two-fields Gaussian Hamiltonian 

H G [p, uj] = ~ (p\w\p) -i(w\p) + ^ ( w I g hs T P~ iuJ o\ \ w ) > ( 63 ) 
and Gaussian averages 

(• • .)g = 1 / VpVw . . . exp (-W G [p, uj]) , (64) 
where the normalization constant A/g is given by 

A/g = / DpDcu exp {-Kg [p, "]) • (65) 

Note that if u> > (attractive case ) eq. (64) makes sense only if the integration 
over the field variable uj is performed in first. With these notations in hands 
we can rewrite eq. (60c) as 

lnS[i/] = lnS MF H + A In A/g + Aln(e* WM ) . (66) 
A cumulant expansion of the last term in the r.h.s. yields the A expansion 

In S = In S MF + A In + A j (6H) G + \ ((SH 2 ) G - (6H) 2 G ) } + 0(A 3 ) , (67) 

where 

[u] = A/g , (68) 

and it will be shown that the third term in the r.h.s. of eq. (67) is of order 
C(A 2 ). 

5.2.1 One-loop approximation 

We want to prove that the CV and KSSHE one-loop approximations for S, 
respectively and EP-\ coincide. Let us first rewrite the formula (56) for 
as the Gaussian average 

E<>] = (exp(i^|Gg)'V)) w . «»> 

Now we focus on [„] = A/q- Performing the (ordinary Gaussian) integra- 
tion over uj first in eq. (65) we find 
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E y > [u] = / Vp exp - Vuj expli (u\p) - - (u\G 



v(2),T 
r HS 



(2), T 
HS 



/Dp 



exp I -- {p\ 



G 



(2),T 
HS 



u>lp 



(70) 



Now we make use of the two following properties of Gaussian integrals 



M a = (71a) 
exp(i^|%>)) =^expQ^|ab>)) , (71b) 



to rewrite 



JX>p exp 



G 



(2),T 
HS 



-1 



/£>p exp ( -\\p\ 



ex P o (P\ W \P) 



G 



(2), T 
HS 



G 



(2), T 
HS 



(72) 



It follows then readily from eqs. (69) and (71b) that = The two 

identities (71a) and (71b) are easy to show in the homogeneous case and, in 
this case, there are an immediate consequence of the properties of the Gaus- 
sian integral (A. 5). In the general, non-homogeneous case, they are derived in 
Appendix B. 



Let us first have a look at the correlations of the field p. 



5.2.2 Two-loop approximation 

At the two-loop level we have to compute the averages (5H [u]) G and (dTi [^] 2 ) ( 
In order to do so we have first to generalize Wick's theorem to this special 
kind of Gaussian average defined in eq. (64). 

Let us first have a look at the correlations of the field p. One has from eq. (64) 



(p(l)...p(n)) ( 



J VpVu p(l) . . . p(n) exp {-H G [p, w]) 
JVpVuj exp (-H G [p, w]) 



JVpp(l) . . .p(n)exp 


(-*(> 




^HS 


-1 

— w 




JVp exp 


hi 


p 


^HS 


1 -1 

— w 
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where we have performed the (ordinary) Gaussian integral on the field u first. 



^(2), Tl - 1 
^HS 



— W 



We are thus led to the ordinary Gaussian average (provided that 
is positive definite, which we assume), therefore the usual Wick's theorem is 
applied which yields 



<P(1)> 



T 



(p(l)p(2))5 



^HS 



— W 



(p(l)...p(n))J = 0(forn>3). 



(73) 



We take now advantage of the general relations (33) between the truncated 
correlations of p and iq in the CV formalism to get 



M1))g = 

(o;(1M2))S = -w(l, 1>(2, 2')G© T (1', 2') = -A, (1, 2) 
(w(l)...w(n))5 = (for n > 3) , (74) 



from which Wick's theorem follows 



(u(l) . ..u{n)) G = in odd ) , 

(w(l)...o;(n)> G = (-l) B £ II A*o(M 2 ) (n even ) , (75) 

pairs 

which was an expected result thanks to the correspondence ipx z^cv- We 
have now at our disposal all the tools to compute 



(SH M> G = ^A y dl . . . d4 Gffi [17 + po] (1, . . . 4) 
x(^(l)...c(4)) G + 0(A 2 ) 

= ^|dl...d4Gffi[l7 + ¥ ,o](l,...4) 

xA W) (l,2)A W) (3,4) + G(A 2 ), (76) 

where we have made use of Wick's theorem (75). We note that (<57i[u;]) G 
will not contribute to InS at the two-loop order and it remains to compute 
{SH 2 [u]) G . One finds 
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(5H 2 [ w ]) g = — - Jdl... d3'G<£l [V + <p ] (1, 2, 3)Gg [* + y? ] (l', 2', 3') 
x(a;(l)...a;(3')} G + 0(A 2 ) 

= \jdl... dS'G^l [V + y? ] (1, 2, 3)G® [V + y? ] (1', 2', 3') 
x{iA V0 (l,2)A V0 (l',2')A V)0 (3,3') 

+ l - A vo (1, 1') A w (2, 2') A vo (3, 3') } + 0(A 2 ) , (77) 

where once again Wick's theorem (75) was used. Gathering the intermediate 
results (76) and (77) one has finally, after inspection 

A {(5H) G + \ (lbU 2 ) G - (5H) 2 G )} = A 2 lnS (2) [u] + G(A 3 ) (78) 

with 

H (2) [„] = [v] . (79) 

We have thus shown that the one and two-loop expressions for In S [u] coincide 
in the KSSHE and CV representations. This coincidence is likely to be true 
at all orders in the loop-expansion. 



6 The pressure and the free energy of the homogeneous fluid 



6. 1 The pressure 



In this section we restrict ourselves to the homogeneous case, therefore In S [u] = 
Vf3P{v), where P denotes the pressure and (3 A [p] = Vf3f(p) where / is the 
Helmoltz free energy per unit volume. The two-loop expression that we de- 
rived for P in Sec. 5 is too much complicated to be of any use in practical 
calculations since it involves the 3 and 4 density correlation functions of the 
HS fluid which are unknown, whereas G^s is known approximatively, for in- 
stance in the Percus-Yevick (PY) approximation [12]. A simple but coherent 
approximations for Gjjs and G^g will be proposed now. 

Recall first that it follows from their definitions [26,27] (see e.g. eqs (15)) that 
the Ghs' T M satisfy to the following relations 

Mn + 1) G ™ T M (1 ' ' ' ' ' ^ = G<SS+1) ' T M (1> ■■■>»>» + 1) ■ ^ 
For a homogeneous system (in which case v is a constant) one infers from this 
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equation that 

d 



d\...dn G^ 1] ' T M (1, • • • , n, n + 1) = j^p m (u) = p®> (y) , (81) 
where pns iy) is the number density of hard spheres at the chemical potential 



v. 



In the rest of the section we will adopt the following approximation 

G$ T [u] (1, . . . , n) = p& 1] (is) 5(n, 1) . . . 5(2, 1) for n > 3 . (82) 

Note that this hypothesis is coherent with the exact relations (80) and (81). 
However, for n = 2, we adopt for G^' T [i/] (1,2) some known approximation 
of the theory of liquids (PY approximation for instance). The free propagator 
has then a non trivial expression. In Fourier space it reads 

**o(g) = 1M , (83) 

l-w(q)G^ [v + (po](q) 

where 

GSl' T (q) = Jd"xe-^G^ T (x) (84) 

( 2'} T 

is the Fourier transform of G H g [is] (x). The set of approximations that we 
have just discussed is reasonable as long as the range of the KSSHE field 
correlation functions is (much) larger than than the range of the HS density 
correlation functions. This would be true if w is the long range pair interaction 
or, in general, in the vicinity of the critical point. 

With the hypothesis (82) it is easy to obtain the two-loop order approximation 
for the pressure. One finds 



(3P(v) = f3P(°\v) + \(3P^(is) + \ 2 pP$\v) + 0(\ 3 ) , 

2 

(3PM(is) = pP MF (is) = P m (V + Vo) - ^ , 
PP {1) (v) = -l I ^(l-w( q )G^ T [v + ipo ]( q )) , 



(3) 

^PP)( 1/ ) = ^_AJ (0) + 



(2) 



+ 



fo 
(2)T2 

Po 



-A W ,(0)AJ ( ) 



12 



d 3 x A 3 (x) , 



(85) 



where = (V + <Po) and / q = / d 3 q/(2nf. 
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6.2 The free energy 



In order to compute the free energy f(p) at the two-loop order we need the 
expression of the density only at the one-loop order (this a well known prop- 
erties of Legendre transform [2] that will be checked explicitely further on). 
We thus define 

p(u) = ^P = p + Ap(u) 
Ap(v)=\pW(v) + 0(\ 2 ), (86) 
where p = p H s 7 + <fo) and 



= f/W0>(l + £), (87) 

where, in order to derive the last line, we made use of eqs. (80) and (82). 
In order to obtain the expression of dipo/du one derives the homogeneous 
stationary condition (38) with respect to v which gives us 



dipo 

dv 1 - p ( o ] w(0) 
Introducing this result in the last line of eq. (87) we find finally that 



= A m 0) • m 

2 l-p^(0) 

We are now in position to compute the free energy 



(3f{p) = pv-(3P{v) 

= pu- (3P {0) (y) - A/3P (1) {v) - \ 2 [3P {2) (u) + C(A 3 ) . (90) 

Our task is now to reexpress the r.h.s. as a function of p = p + Ap(u) which 
will be done along the same lines that those used in [18]. We shall compute 
separately three contributions to f3f. 

i) Computation of X = pv — (3P(°\v) 

We first note that X can be rewritten as 

X = Po u - pp(°\is) + vAp = (3f {0 \p ) + vAp , (91) 
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where / (0) (Po) = /mf(po) is the MF free energy of the homogeneous system 
at the density p , i- e - (cf. eq. (42)) 

Pf {0) (po) = Pfm(po) - lp w 2 (0) + ±p w(0) . (92) 
In order to relate f^(po) to f(p] we perform a second order Taylor expansion 



Pf (0) (p = Po + Ap) = Pf^(po) + vAp - \c MF (0) [Ap] 2 + G(A 3 ) , (93) 

where the Fourier transform at q = of the MF two-body vertex function is 
computed from eq. (44b) with the result 

CW(0) = -1/G HS (0) + iS(0) = . (94) 

Po 

Combining the intermediate results (92), (93), and (94) one has finally 



X = [3f^(p) + ±C MF (0) [Ap] 2 + 0(X 3 ) 



2 



= (3f (0) (p) ~ 4-A°(0) nw /" kv—, + 0{\ 3 ) . (95) 



Po 



(2) 



2 



8 ^ 7 P «(l-p^(0)) 



We note that, as claimed at the beginning of the paragraph, the knowledge of 
Ap at the one-loop order is sufficient for our purpose. 

ii) Computation of pP^\v) 

As can be inferred from eq. (85) (3P^ l \v) may be seen as a function of the sole 
variable p = Phs(v+<Po) since for the HS reference system there is a one to one 
correspondence between densities and chemical potentials -at least in the fluid 
phase, away from the liquid/solid transition-. Therefore (3P^\v) = pP^(po), 
with some abusive notations. One thus have 



PPM( P = po + Ap) =/3PV(po) + ^^M Ap + o(X) , 

opo ov 

= ^ (1) (Po) + |^^ + 0(A), (96) 
opo A 

where the last line follows from eqs. (86) and (87). At this point we notice 
that 

dpo ^ T ( m Po 1} i / dv 

* = Gmf( ' = 0) = = 11 W» • (97) 
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which allows us to write finally 



)(!)/ 



>(i). 





r (2)" 

Po 


2 


(1 


- Po^(0)) 



(3P (1) (p) = -lJln(l-w( q )G^ p ( q )) , 



(98a) 
(98b) 



where, in the second equality, the subscript emphasizes that the truncated HS 
pair correlation function has to be computed at the density p. 



The computation of f(p) is now completed and gathering the intermediate 
results one has finally, setting A = 1 



+ 



\f {ln(l-w( q )Gl^ p ( q ))+pw(q)} 



(3) -I 

^AJ(O) + U%0) 



(2) 

Po 



(2)1 2 

Po 



n {1) 

Po 



12 



J d 3 x A 3 p {x) . (99) 



Some comments on this result are in order. 



i) The last line in eq. (99) gathers the two-loop contributions of (3f(p). At this 
order of the loop expansion it is thus legitimate to replace A Vo by A p where 
the propagator is evaluated at the density p rather than p since p—po = 0(A). 
More explicitly 

w(q) 



A„(?) = 



1 



Ggs,„(?)tf(?) 



(100) 



ii) Similarly the derivatives which enter the last line of eq. (99) can be 
evaluated at the density p rather than po- One can thus write, taking advan- 
tage of the one to one correspondence between the HS densities and chemical 
potentials 



Po 1} = 



^ (P) 



(2) 

Po 



dpi 



(i) 



dv 



-v, 



(2) 
HS 



(P) 



4kp) 



l 3 



(3) 

Po = 



dp'? 


3 


~4Up)_ 


2 - 4 3 S(P>HS(P) 


dv 






~4kp) 


5 



(101) 
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where ^hs (p) denotes the nth derivative of the HS chemical potential with 
respect to the density (it can be computed in the framework of the PY or 
Carnahan-Starling approximations for instance [12]). 

iii) It must be pointed out that, quite unexpectedly, the reducible diagram 
has not be cancelled by the Legendre transform. Usually, in statistical field 
theory it is the case (cf. [2]). The reason is that the chemical potential v is 
not the field conjugate to the order parameter m = (<p) K of the KSSHE field 
theory. However one of us have shown elsewhere [18,23] that for the symmetric 
mixtures of charged hard spheres only irreducible diagrams contribute to (3f^ 2 \ 

iv) All the quantities which enter eq. (99) can be computed numerically (for 
instance in the PY approximation). We plan to test the validity of our ap- 
proximation in future work. 



7 Conclusion 

Using the CV method we reconsider the basic relations of statistical field 
theory of simple fluids that follow from this approach. In contrary to the 
KSSHE theory [11] the corresponding CV action depends on two scalar fields 
- field p connected to the number density of particles and field uo conjugate to 
p. Explicit expressions that allow to relate between themselves the correlation 
functions defined in different versions of the theory are derived. 

For a one-component continuous model of fluids, consisting of hard spheres 
interacting through the short-range pair potential, we calculated the grand 
partition function in a systematic way for both versions of statistical field 
theory using the loop expansion technique. As it was expected, in all the 
orders of loop expansion considered, both versions of the theory produced 
indeed the same analytical results. The expressions for the pressure as well as 
for the free energy are derived in a two-loop approximation that is the next 
step in comparison with the results obtained recently by one of us [11] within 
the KSSHE theory. In fact this is a new type of approximation which we plan 
to test in our future work. 

In contrast to the usual statistical field theory [2], our results demonstrate 
that for the case of simple fluids the reducible diagram is not cancelled by 
the Legendre transform. It is due to the fact that our field theory is a non- 
standard in the sense that the coupling between the internal field ip and the 
external field v is non-linear yielding new expressions for (3f(p). 

From our analysis of the CV and KSSHE transformations we can also con- 
clude that the former has some important advantages which could be very 
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useful for more complicate models of fluids. In particular, it is valid for an 
arbitrary pair potential (including a pair interaction w(l,2) which does not 
possess an inverse) and is easily generalized for the case of n-body interparticle 
interactions with n > 2. 
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A Functional measures and integrals 



In this appendix we give some details on functional measures and integrals. Let 
us consider a real scalar field y(x) defined in a cube C3 of side L and volume 
V = L 3 . We assume periodic boundary conditions, i.e. we restrict ourselves to 
fields which can be written as a Fourier series, 

where A = (2tt/L) Z 3 is the reciprocal cubic lattice. The reality of tp im- 
plies that, for q 7^ </5_ q = </?*, where the star means complex conjugation. 
Following Wegner [19] we define the normalized functional measure Vip as 



difqdip-q = 2 d^lipq d^sif^ for q 7^ . (A. 2b) 

Eq. (A. 2) can be conveniently rewritten as 

Vip = n M ^ ^ , (A.3) 

where the sum in the r.h.s runs over only the half A* of all the vectors of the 
reciprocal lattice A (for instance those with q x > 0). With the definition (A. 2) 
we have 
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AC, = 



J £V exp ^~ (y\w 1 \ip 



= ex P H ln M<?) 



L— >oo 
> i 



exp J ]nw(q)^ , (A.4) 



where w(l,2) is positive and satisfies w(l,2) = w(2, 1) = w(x\2). It is worth 
noting that in this case we have trivially AC- 1 = 1/Af w . More generally and 
with the same hypothesis we have the useful identity 

(exp ((i) ((f\uj))) w =Af- 1 J V<p exp (-± (flw' 1 ^) + (i) H^)) 

= exp(+(-)^(wHo;)j , (A.5) 
where a; is a real scalar field. 

We define now the "functional delta" distribution 5jr [A] as [1] 

5r [A] = Jvuj exp(i(w|A)) , (A.6) 
where both cj and A are rea/ scalar fields defined on C. Since 

HA) = / rf 3 x c(x)A(x) = "-q*q , (A.7) 

it follows from eq. (A. 2) that we have more explicitly 

b T [A] = V2^V 5 (Ao) [] [* v 6 ( KA q) 6 (^q)] ■ ( A - 8 ) 

qeA* 

Therefore 

J V\ 5 T [A] = I , (A.9) 

and, more generally 

\v\T [A] Sjr [A - Ao] = T [Ao] , (A. 10) 
where T [A] is some arbitrary functional of the field X(x). 

B Two useful identities 



We prove here two identities used in Sec. 5. Firstly let if (x, y) be a definite 
positive operator assumed to be symmetric. The variables x and y take their 
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values in the cube C3 with periodic boundary conditions. We define (provided 
that it exists) 

M K = J Pep exp (-1 , (B.l) 

and we want to prove that M K -i = 1 /Mr- If K(~k, y) = K(x — y) = if (y — x), 
i.e. in the homogeneous case, it is trivial as pointed out in Appendix A. In the 
general case we still have [2] 

J Vtp exp (-± (iflK-^ip) + i (J|^>) = M K exp (~\ (J\K\J)} , (B.2) 
where J is some real scalar field. Therefore 



M K -i = J VJ exp (-±(J\K\J?) 

= JuS VJVv exp i~l i^' 1 ^) + 1 

= J^JVp exp (-± (^|if-V)) ^ M 



Mr ' 

where we have made use of the representation (A. 6) of 8? [ip]. 



(B.3) 



Secondly let us consider two definite positive symmetric operators a(x, y) and 
6(x, y). We define (provided that these objects exist) 

(I, , . ,\\ /^exp(-i(^|6-» + |(^|a|^)) 
exp -^ab ) = — ^ — — — — , ^, (B.4) 

and we want to prove that 

exp (± (cp\a\^ ) 6 = (exp Q ^ . (B.5) 

Note that these expressions are well defined only if (a— b~ r ) > and (b— or 1 ) > 
0, i.e. afr < 1. If a(6)(x, y) = a(6)(x — y) = a(6)(y — x) (i.e. in the homogeneous 
case) the result is a trivial consequence of eq. (A. 5). 

In the general case one proceeds as follows. Since a > and b > then one 
can define the operators a ±l l 2 and b ±l l 2 which are all symmetric and positive 
definite. Let us perform the change of variable 

<p(l)=W(l,l')a- 1 '*(l'XMl") (B.6) 

in eq. (B.4). The functional Jacobian of the transform J = det 5<p(l) / 5if}(2) = 
det 6 1 / 2 (1, l')a -1 / 2 (l', 2) is a non zero constant which can be picked out of the 
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integral. Moreover it is easy to show that 



(^a" 1 ]^), (<p\a\<p) = (V>|£#> 



(B.7) 



therefore one has 



( 




Jib 



J/^exp(-l^|a^|^) + |(^|^)) 
J/^exp(-i(^|a-^)) 




(B.8) 
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